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SUMMARY/OVERVIEW 

The  objeetive  of  the  proposed  researeh  is  to  understand  the  roles  of  loeal  limit  phenomena,  e.g. 
loeal  extinetion  and  re-ignition,  interaetions  between  flow  eompression  and  fast-reaetion  indueed 
dilatation  (reaetion  eompression),  and  to  identify  the  rate-limiting  fuel  oxidation  pathways  in 
high-speed  turbulent  flames.  To  aehieve  this  objeetive,  systematie  eomputational  diagnosties 
were  developed  and  performed  based  on  elementary  flames  and  turbulent  flames  simulated  with 
direet  numerieal  simulations  (DNS).  During  the  eurrent  reporting  period,  progress  has  been  made 
in  the  following  tasks:  1)  Computational  diagnosties  were  developed  to  identify  loeal  ignition 
and  extinetion  and  reaetion  fronts  in  turbulent  flames  simulated  using  DNS.  2)  Effeets  of 
surrogate  jet  fuel  eomposition  and  fuel  eraeking  reaetions  on  high  temperature  eombustion  were 
studied  using  a  bifureation  analysis.  3)  Redueed  meehanisms  for  high-temperature  eombustion  of 
n-butane  and  n-dodeeane  with  lumped  fuel  eraeking  reaetions  were  developed  and  validated.  4) 
2-D  DNS  was  performed  to  investigate  the  fuel  eraeking  effeets  in  high  temperature  and  strongly 
turbulent  premixed  flames.  5)  A  highly  effieient  ODE  solver,  dynamie  adaptive  hybrid 
integration  (AHl),  was  developed  for  stiff  ehemistry. 
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TECHNICAL  DISCUSSION 


I.  Chemical  explosive  mode  analysis  (CEMA)  for  computational  flame  diagnostics 

The  method  of  chemical  explosive  mode  analysis  (CEMA)  is  a  systematic  approach  to 
identify  limit  flame  phenomena,  including  local  ignition,  extinction,  and  premixed  and  non- 
premixed  reaction  fronts,  in  both  laminar  and  turbulent  flames.  A  chemical  explosive  mode 
(CEM)  is  an  eigenmode  of  the  chemical  Jacobian  associated  with  a  positive  eigenvalue.  CEM  is 
a  chemical  property  that  distinguishes  a  combustible  mixture  from  other  substances,  and  the 
interaction  between  CEM  and  the  mixing  processes  can  result  in  local  ignition  and  extinction. 
Therefore  CEMA  is  an  effective  diagnostics  to  identify  local  limit  phenomena. 

In  the  present  project,  CEMA  is  used  to  identify  and  track  the  local  ignition  and 
extinction  in  turbulent  flames  simulated  by  DNS  at  Sandia.  Two  DNS  datasets  previously 
generated  at  Sandia  are  employed  for  the  analysis:  a  non-premixed  temporal  jet  flame  of 
ethylene-air  diluted  with  N2  (DNS  by  D.O.  Lignell),  and  a  premixed  turbulent  counterflow  flame 
of  hydrogen- air  supported  by  super- adiabatic  hot  products  (DNS  by  S.  Lyra).  The  two  DNS 
datasets  show  significantly  different  flame  extinction  behaviors:  the  non-premixed  ethylene  jet 
flame  is  characterized  by  a  folded  S-curve,  which  is  typical  in  practical  combustors,  and  the 
counterflow  flame  features  a  stretched  S-curve  with  the  flame  being  support  by  adiabatic  hot- 
product  (5 OK  above  adiabatic  flame  temperature).  Both  the  DNS  datasets  were  generated  from 
Sandia’ s  multi-million-CPU-hour  supercomputing  and  are  high  fidelity  data  sources  for 
computational  diagnostic  benchmarking  and  turbulent  combustion  model  creation  and  validation. 

The  CEMA  result  for  the  non-premixed  temporal  jet  flame  for  ethylene/air  is  shown  in 
Eig.  1  for  the  center-cut  plane  of  the  3-D  DNS  at  t  =  0.1ms.  Extensive  local  extinction  can  be 
observed  on  the  temperature  plot  as  indicated  by  the  discontinuities  on  the  non-premixed  reaction 
fronts  (high-T  layer  indicated  by  the  dark  red  colors).  On  the  CEMA  plots,  the  near-  and  post- 
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extinction  mixtures  are  unambiguously  captured  as  the  red  spots  (with  CEM)  on  the  eigenvalue 
plot  and  the  yellow-red  spots  on  the  Da  plot,  as  mixtures  becomes  explosive  at  near-extinction 
conditions  where  the  incomplete  combustion  results  in  high  concentrations  of  reactants  while  the 
temperature  is  still  high.  Compared  with  the  temperature  based  criteria,  the  CEMA  based 
approach  is  more  robust  and  universal  to  capture  the  local  extinction  (and  possible  re-ignition). 
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Eigure  1:  A  2-D  snapshot  of  a  3-D  DNS  of  a  non-premixed  temporal  jet  of  N2-dilited  ethylene-air 
experiencing  substantial  local  extinction.  The  local  extinction  was  captured  by  CEMA  on  the 
eigenvalue  plot  (as  red  spots)  and  the  Da  plot  (yellow-red  spots). 
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Figure  2:  A  2-D  snapshot  from  a  3-D  DNS  of  a  premixed  eounterflow  flame  for  hydrogen/air 
supported  by  super-adiabatie  hot  produet  (red  zones  on  the  temperature  plot).  A  flame  hole  is 
visualized  by  CEMA  on  the  eigenvalue  plot  (with  the  dark-blue  post-reaction  zone  missing)  and 
the  species  explosion  index  (El)  plot.  The  controlling  reactions  are  further  identified  by  the 
reaction’s  participation  index  (PI). 


The  CEMA  results  on  the  eounterflow  flame  of  premixed  hydrogen/air  are  shown  in  Eig. 
2.  It  is  noted  that  the  hydrogen  flame  is  supported  by  super-adiabatic  hot  products  (5 OK  above  the 
adiabatic  flame  temperature),  such  that  local  flame  extinction  is  not  discemable  from  the 
temperature  plot.  In  contrast,  CEMA  successfully  identified  the  flame  hole  (local  extinction) 
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through  the  identification  of  the  missing  post-fiame-zone  in  the  flame  hole  as  shown  on  the 
eigenvalue  plot  of  Fig.  2.  It  is  seen  that  the  self-sustaining  premixed  flame  section  features  a 
strong  post-reaction  zone  structure  as  shown  in  dark  blue  colors,  while  the  flame  hole’s  post¬ 
reaction  zone  is  mostly  suppressed  (thus  the  flame  is  not  self-sustaining)  as  indicated  by  the  light 
blue  colors.  The  flame  hole  is  further  evident  on  the  species  explosion  index  (El)  plot.  The  El 
value  of  a  species  indicates  the  normalized  contribution  of  a  species  to  the  local  explosive  mode. 
It  is  seen  that  the  controlling  species  (those  with  large  El  values)  are  dramatically  different  for 
strongly  burning  flames  vs.  the  flame  hole,  showing  that  the  controlling  chemistry  will  change  in 
a  flame  once  local  extinction  occurs. 

The  CEMA-based  criteria  are  therefore  shown  to  be  effective  and  can  be  extended  as 
computational  diagnostics  to  systematically  detect  and  track  local  ignition,  extinction  and 
reaction  fronts  in  complex  turbulent  flame  fields,  providing  a  robust  local  flame  marker  for 
switching  between  different  turbulent  combustion  models  that  are  suitable  for  different  local 
flame  features. 

2.  Effects  of  surrogate  components  on  jet  fuel  ignition  and  extinction  and  a  bifurcation  analysis 

Jet  fuels  involve  a  large  number  of  components  with  different  molecular  structures,  e.g. 
alkanes,  alkenes  and  aromatics,  primarily  ranging  from  C7  to  Cl 6.  It  is  computationally 
expensive  to  solve  such  a  large  number  of  fuel  species  and  the  myriad  intermediate  species  and 
reaction  pathways  associated  with  these  fuel  components  in  turbulent  flame  simulations.  As  such, 
surrogate  mixtures  with  one  or  a  few  surrogate  components  are  frequently  employed  to  develop 
chemical  kinetic  models  to  capture  the  major  properties  of  the  real  fuels,  e.g.  the  distillation 
curves,  thermodynamic  properties,  ignition  and  extinction  behaviors,  and  flame  speeds.  In  the 
present  project,  the  viability  of  using  a  single-component  surrogate  to  capture  high- temperature 
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jet  fuel  combustion  features  (ignition  and  extinction)  was  computationally  explored,  with  the 
results  explained  using  a  rigorous  bifurcation  analysis. 

Three  types  of  hydrocarbons,  including  normal-alkanes,  cyclo-alkanes,  and  aromatics  are 
considered  as  primary  jet  fuel  surrogate  components  in  the  present  study.  It  was  found  that 
normal-alkanes  with  different  chain  lengths  (C8~C12)  and  branched  cyclo-alkane  (represented  by 
n-butyl  cyclohexane)  show  similar  ignition  and  extinction  behaviors  in  high-temperature 
combustion  (e.g.  above  lOOOK).  As  such  pure  n-dodecane  can  be  an  adequate  representative 
component  for  the  normal-  and  cyclo-alkanes  components  in  jet  fuels  in  terms  of  capturing  high- 
temperature  ignition  and  extinction  behaviors. 


1000/T,  1/K 


Figure  3:  (a)  Ignition  delay  time  as  a  function  of  initial  temperature  in  constant -pressure  auto¬ 
ignition,  and  (b)  the  ^-curves  of  perfectly  stirred  reactors  (PSRs),  for  n-dodecane,  toluene,  and 
binary  mixtures  of  n-dodecane  and  toluene  with  0,  20%,  and  50%  toluene  in  mole,  respectively. 


It  was  further  shown  in  Fig.  3  that,  while  pure  aromatics  (e.g.  toluene)  features 
significantly  ignition  and  extinction  properties  from  n-dodecane,  blending  a  small  to  moderate 
amount  of  toluene  (say  up  to  50%  in  mole)  in  n-dodecane  only  induces  minor  changes  in  the 
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ignition  and  extinction  states.  This  result  indicates  that  pure  n-dodecane  can  be  a  sufficiently 
accurate  jet  fuel  surrogate  for  high-T  ignition  and  extinction,  such  that  a  simple  chemical  kinetic 
mechanism  for  pure  n-dodecane  can  be  adequate  to  accurately  predict  high-T  ignition  and 
extinction  of  jet  fuels. 

To  explain  the  underlying  reasons  for  this  observation,  a  bifurcation  analysis  is  employed 
to  identify  the  controlling  species  and  reactions  for  jet  fuel  ignition  and  extinction.  The 
bifurcation  analysis  is  focused  on  (quasi-)  steady  flames  that  feature  the  combustion  S -curves. 
The  upper  turning  point  of  the  S-curve  is  corresponding  to  the  extinction  state  and  the  lower 
turning  point  to  the  ignition  state.  The  extinction  and  ignition  states  are  important  for  flame 
stabilization  and  are  mathematically  described  as  bifurcation  points  associated  with  zero-crossing 
eigenvalues  of  the  Jacobian  matrix.  The  controlling  reactions  of  the  ignition  and  extinction  can 
therefore  be  identified  by  quantifying  the  contribution  of  each  species  and  reaction  to  the  zero¬ 
crossing  of  the  eigenvalue,  and  such  normalized  contribution  from  a  reaction  to  the  bifurcation  is 
defined  as  the  bifurcation  index  (BI),  which  is  used  in  the  present  study  to  subsequently  identify 
the  controlling  reactions  for  ignition  and  extinction. 

The  reactions  and  the  mixing  process  with  large  BI  values  are  listed  in  Figs.  4a  and  b  for 
the  extinction  and  ignition  states,  respectively.  Fig.  4a  shows  that  the  extinction  states  of  the 
different  binary  mixtures  with  different  amounts  of  toluene  addition  are  controlled  by  the  mixing 
process  and  a  set  of  reactions  involving  small  molecules  that  are  relevant  to  CO  formation  and 
the  hydrogen  chemistry,  while  the  reactions  for  fuel  cracking  for  either  n-dodecane  or  toluene  do 
not  play  a  significant  role  at  the  extinction  state.  At  the  ignition  state,  it  is  shown  in  Fig.  4b  that 
the  controlling  reactions  for  the  different  mixtures  primarily  involve  ethylene,  vinyl  radical,  alkyl 
radical  and  other  small  species  that  are  important  intermediate  products  of  beta-scission.  The 
results  in  Fig.  4  indicate  that,  since  fuel  cracking  through  beta-scission  becomes  exceedingly  fast 
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at  high  temperatures,  they  are  not  rate-limiting  at  either  the  extinction  or  the  ignition  states.  As 
such,  it  is  possible  to  capture  high-T  jet  fuel  ignition  an  extinction  behaviors  using  lumped  fuel 
cracking  reactions,  such  that  the  jet  fuel  chemistry  can  be  significantly  simplified  and  expedited 
for  time  integration  at  high-T  flame  simulations. 


H  n-dodecane 

H  80%  n-dodecane+20%  toluene 
H  50%  n-dodecane+50%  toluene 


HC0+02  =  CO+H02 

a  C6H5CH3  +  02  =  C6H5CH2  +  H02  ,,,,,, jj;! 

CH3+H(+M)  =  CH4{+M) 

?  C6H5CH2  +  H02  =  C6H5CHO  +  H  +  OH  bgp, . 

HCO+H  =  CO+H2  ^ 

H02+0H  =  H20+02 

H+02(+M)  =  H02{+M)  g 

CH3+H02  =  CH30+0H 

H+OH+M  =  H20+M  ™ 

OH+OH(+M)  =  H202{+M) 

H+02  =  O+OH 

aC3H5+H02  =  0H+C2H3+CH20 

CO+OH  =  C02+H 

C2H4+OH  =  C2H3+H20 

Mixing 

Mixing 

(a)  Bifurcation  index  at  extinction  (b)  Bifurcation  index  at  ignition 

Figure  4;  Reactions  ranked  in  bifurcation  index  (BI)  for  the  (a)  extinction  and  (b)  ignition  states 
on  the  5'-curves  of  perfectly  stirred  reactors  (PSRs)  in  Fig.  3b  for  the  binary  mixtures  of  n- 
dodecane  and  toluene. 


3.  Reduced  mechanisms  for  n-butane  and  n-dodecane  with  lumped  fuel  cracking  reactions 

Detailed  chemistry  is  an  integral  component  of  turbulent  flame  simulations.  As  such  an 
accurate  and  efficient  chemistry  description  is  critical  for  CFD  predictions.  To  develop  highly 
efficient  yet  still  accurate  reduced  mechanisms  for  jet  fuel  combustion,  the  above  observation  for 
the  fuel  cracking  is  exploited  and  the  fuel  cracking  reactions  are  lumped  (in  collaboration  with 
Hai  Wang)  to  significantly  reduce  the  amount  of  intermediate  species  generated  during  fuel 
cracking.  Two  fuels,  namely  n-butane  and  n-dodecane,  are  selected  for  the  reduction:  n-butane 
was  selected  as  it  is  the  smallest  alkane  that  features  clear  fuel  cracking  pathways,  while  n- 
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dodecane  was  selected  as  it  is  an  important  jet  fuel  surrogate  for  high-T  flame  conditions  as 


discussed  above. 


n-butane/air 


1000/(T,  K) 


Figure  5;  Validation  of  the  lumped  and  reduced  mechanisms  for  n-butane/air:  a)  laminar  flame 
speed,  and  b)  ignition  delay  time  at  different  pressures  and  equivalence  ratios.  Solid  lines: 
detailed  USC-Mech  II;  dashed  lines:  lumped  fuel  cracking  reactions;  symbols:  28-species 
reduced. 

The  detailed  USC-Mech  II  mechanism  with  III  species  is  used  as  the  starting  mechanism 
for  n-butane.  The  fuel  cracking  reactions  were  first  lumped  (by  Hai  Wang  and  student),  resulting 
in  a  58-species  mechanism.  Systematic  mechanism  reduction  methods  based  on  directed  relation 
graph  (DRG),  sensitivity  analysis  and  quasi  steady  state  approximations  (QSSA)  were  further 
employed  to  reduce  the  mechanism  to  28  species.  The  detailed,  lumped  and  reduced  mechanisms 
are  compared  in  Fig.  5  for  flame  speed  and  ignition  delays,  and  close  agreement  was  observed  for 
the  displayed  and  other  flame  conditions. 

For  n-dodecane,  a  JetSurF  mechanism  (from  Hai  Wang’s  group)  with  123  species  and 
lumped  fuel  cracking  reactions  was  employed  as  the  starting  mechanism.  The  reduction  covered 
auto-ignition  and  perfectly  stirred  reactors  for  equivalence  ratio  range  of  0. 5-1.5,  initial 
temperature  higher  than  1000  K,  and  at  atmospheric  pressure.  The  reduced  mechanism  consists 
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of  only  24  species,  with  the  high  accuracy  of  the  reduced  mechanism  shown  in  Fig.  6.  The 
effectiveness  to  use  lumped  fuel  cracking  steps  to  simplify  high-T  jet  fuel  chemistry  is  therefore 
evident. 


Residence  time,  s 

10"*  10^"  10^"  10'  10° 


Figure  6;  Ignition  delay  as  a  function  of  initial  temperature  in  constant-pressure  auto-ignition 
(bottom-left)  and  temperature  as  a  function  of  residence  time  in  PSR  with  inlet  temperature  of 
300  K  (top-right),  for  stoichiometric  n-dodecane/air  at  atmospheric  pressure,  calculated  with  the 
detailed  and  reduced  mechanisms,  respectively. 


4.  2-D  DNS  of  premixed  turbulent  n-butane  flames 

To  further  study  the  effects  of  fast  turbulent  mixing  on  fuel  cracking  reactions,  2-D  DNS 
for  premixed  n-butane/air  flames  (equivalence  ratio  =  0.6)  were  performed  in  the  thin  reaction 
zones  regime  and  the  broken  reaction  zones  regime,  respectively.  The  flames  are  initialized  with 
the  1-D  premixed  flame  profile  and  isotropic  turbulence,  and  the  flames  propagate  into  fresh 
mixtures  at  500K  at  a  pressure  of  5  atm,  with  the  flame  fronts  disturbed  by  the  turbulence 
fluctuations.  The  computational  domain  is  10mm  by  15mm  with  5  uniform  grid  points. 
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Figure  7;  DNS  of  a  2-D  premixed  n-butane/air  flame  in  the  broken  reaction  zones  regime:  a)  the 
fuel  concentration  and  b)  the  concentrations  of  selected  fuel  cracking  products  at  t  =  0.7ms.  The 
white  isolines  are  the  premixed  reaction  fronts  (eigenvalue  zero-crossing)  identified  by  CEMA. 

The  flame  structure  at  t  =  0.7ms  for  the  case  in  the  broken-reaction  zones  regime  is 
shown  in  Fig.  7.  The  DNS  results  show  that  despite  the  strong  turbulent  mixing,  fuel  molecules 
have  never  been  observed  to  penetrate  into  the  reaction  zone  of  the  flames  indicated  by  the  white 
isolines  in  Fig.  7,  and  the  fast  fuel  cracking  assumption  is  valid  for  flames  even  with  strong 
turbulent  mixing.  Further  statistical  analysis  of  the  DNS  data  shows  that  the  fuel  cracking 
reactions  are  mostly  complete  at  temperature  between  1500  K  and  1600  K,  as  such  the  lumped 
fuel  cracking  reactions  do  not  have  significant  effects  on  the  chemistry  in  the  reaction  zones  in 
the  strongly  turbulent  premixed  flames. 

5.  Fuel  cracking  in  near-extinction  flames 

To  further  study  the  validity  of  fuel  cracking  lumping  under  extreme  combustion 
conditions,  such  as  in  near-extinction  flames,  the  flame  responses  to  the  flow  strain  was  tested  in 
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1-D  counterflow  non-premixed  and  premixed  flames  using  the  detailed  JetSurF  1.0,  lumped- 


detailed  model  (JetSurF  1.0-1)  and  the  24-speeies  redueed  n-dodecane  model. 
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Figure  8;  (a)  Comparison  of  the  maximum  temperature  Fmax  in  counterflow  non-premixed  flames 


as  a  funetion  of  the  strain  rate,  calculated  with  the  detailed,  lumped-detailed  and  24-species 


reduced  models;  (b)  Normalized  total  mass  fraetion  of  C4+  speeies,  and  that  of  the  primary  fuel 


cracking  products,  including  CFI4,  C2FI4  and  CsHg,  at  three  different  strain  rates  as  predicted  by 


the  detailed  model.  The  flames  are  at  atmospheric  pressure  and  an  inlet  temperature  of  300  K 
with  the  fuel  jet  comprised  of  50%  (mol)  n-dodeeane  in  N2  and  air  at  300  K  as  the  oxidizer  jet. 


Figure  8a  shows  the  maximum  temperature,  eomputed  for  the  non-premixed  counterflow 
flames  as  a  funetion  of  the  strain  rate  02.  It  is  seen  that,  the  lumped  model  agrees  well  with  the 
detailed  model  over  the  entire  curves,  including  the  turning  point,  which  corresponds  to  the 
extinction  state  of  the  flames.  In  worst  ease,  the  24-species  reduced  model  gives  an  error  about 
20%  in  the  extinction  strain  rate.  Figure  8b  shows  the  normalized  total  mass  fraction  of  species 
with  four  or  more  earbon  atoms,  C4-f  and  the  primary  fuel  eraeking  products,  CH4,  C2H4,  and 
C3H6.  The  three  profiles  shown  in  Figure  8b  were  ealeulated  with  rather  different  strain  rates, 
spanning  from  strongly  burning  to  extinction  states.  The  fuel  eraeking  behaviors  are  quite  similar 
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for  all  these  strain  rates.  Specially,  the  C4+  species  primarily  vanish  by  about  1500  K,  indicating 
that  the  fuel  cracking  process  is  completed  in  the  preheat  zone  before  the  main  flame  zone. 
Accordingly,  the  fuel  cracking  products  are  oxidized  in  the  main  flame  zone. 


Figure  9:  (a)  Comparison  of  maximum  temperature  Fmax  in  counterflow  premixed  flames  as  a 
function  of  the  strain  rate,  calculated  with  the  detailed,  lumped-detailed  and  24-species  reduced 
models;  (b)  Normalized  total  mass  fraction  of  C4+  species,  and  that  of  the  primary  fuel  cracking 
products,  including  CFI4,  C2H4  and  CsHg,  at  three  different  strain  rates  as  predicted  by  the 
detailed  model.  The  flames  are  at  atmospheric  pressure  and  an  inlet  temperature  of  300  K  with 
the  twin  jets  composed  of  an  n-dodecane/air  mixture  at  equivalence  ratio  0.7. 

Figure  9  shows  extinction  of  premixed  counterflow  flames  predicted  using  both  the 
lumped-detailed  and  the  24-species  reduced  models.  The  results  are  close  to  that  of  the  detailed 
model,  including  the  near  extinction  states.  Comparison  of  the  premixed  and  non-premixed 
counterflow  results  show  similar  features.  In  all  cases,  fuel  cracking  is  completed  before  1500  K. 
Overall,  the  lumping  approach  is  shown  to  be  valid  for  premixed  and  non-premixed  flames  at 
near-extinction  conditions. 
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6.  A  dynamic  adaptive  hybrid  method  for  stiff  chemistry  integration  (AHI) 

Detailed  ehemical  kinetie  models  are  expensive  to  be  aceommodated  in  CFD  simulations 
due  to  not  only  the  large  sizes  but  also  severe  ehemical  stiffness,  which  renders  the  low  cost 
explicit  time  integration  solvers  infeasible  unless  extremely  small  integration  time  steps  are 
assumed.  As  such,  implicit  solvers  are  typically  employed  in  LES  and  RANS,  and  DNS  of 
incompressible  flows.  However,  fully  implicit  integration  solvers  are  computationally  prohibitive 
for  multi-dimensional  flow  simulations,  and  operator  splitting  schemes  are  widely  used  to 
separate  the  integration  of  chemical  and  flow  processes  into  sequential  sub-steps.  It  was  found 
that,  for  systems  with  significant  radical  sources  from  transport,  splitting  chemistry  from 
transport  may  induce  0(1)  errors.  Similar  splitting  errors  were  also  identified  for  extinction 
problems.  To  solve  the  problem,  a  dynamic  adaptive  hybrid  integration  (AHI)  method  for  stiff 
chemistry  integration  was  developed.  In  AHI,  chemistry  and  transport  are  integrated  together  to 
avoid  the  splitting  errors.  Only  fast  species  and  fast  reactions  are  solved  implicitly  while  the  other 
variables  and  source  terms  are  treated  explicitly  to  reduce  the  computational  cost.  Furthermore, 
an  improved  second  order  AHI  scheme  (AHI2)  was  developed  to  achieve  higher  accuracy. 
Figure  10  shows  the  comparison  of  solutions  of  a  toy  problem  using  AHI  and  the  Strang  splitting 
scheme,  respectively.  The  toy  problem  involves  rate  controlling  reactions  with  0(1)  timescales 
and  a  radical  with  10'^  timescale.  Figure  11  compares  errors  in  the  Strang  splitting,  AHI  and 
AHI2,  respectively,  for  the  toy  problem.  It  is  seen  that  the  Strang  splitting  scheme  gives  0(1) 
errors  unless  the  integration  time  step  is  close  to  or  shorter  than  that  required  for  explicit  solvers, 
while  AHI  accurately  predicts  the  species  concentrations  with  significantly  higher  accuracy. 
AHI2  achieved  even  higher  accuracy  over  a  wide  range  of  integration  time  steps. 
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Figure  10:  Profiles  of  speeies  eonoentrations  in  a  toy  problem,  ealeulated  with  the  first-order  AHI 
method  (open  symbols)  and  the  Strang  splitting  seheme,  respeetively,  in  eomparison  with  the 
exaet  solution  (lines).  Time  steps  for  Strang  splitting  and  AHI  are  both  10'^. 
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Figure  1 1 :  Error  eontrol  behavior  for  Strang  splitting,  first-order  AHI  and  seeond-order  AHI  in 
the  toy  problem.  The  dotted  line  is  a  trend  line  with  a  slope  of  2,  indicating  second-order 
accuracy.  The  dashed  line  is  a  trend  line  with  a  slope  of  1 ,  indicating  first-order  accuracy. 
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Figure  12:  Average  CPU  time  per  step  for  fully  implieit  DASAC,  one-stage  explieit  solver,  AHI 
with  dense  LU  deeomposition,  and  AHI  with  sparse  LU  deeomposition,  respeetively. 

In  term  of  eomputational  effieieney,  Fig.  12  further  eompares  the  fully  implieit  DASAC 
solver,  one-stage  fully  explieit  solver,  AHI  with  dense  LU  deeomposition,  and  AHI  with  sparse 
LU  faetorization.  It  is  seen  that  AHI  with  dense  LU  ean  aehieve  0(10)-0(100)  speedup  faetors 
for  meehanisms  of  different  sizes,  while  AHI  with  sparse  LU  faetorization  ean  aehieve 
eomparable  eost  linearly  proportional  to  the  meehanism  size,  being  eomparable  to  the  highly 
effieient  single-stage  explieit  solver. 
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